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ABSTRACT 

The starburst galaxy M82 contains two ultraluminous X-ray sources (ULXs), 
CXOM82 J095550.2+694047 (=X41.4+6Q) and CXOM82 J095551. 1+694045 (=X42.3+59), which 
are unresolved by XMM-Newton. We revisited the two XMM-Newton observations of M82 and 
analyzed the surface brightness profiles using the known Chandra source positions. We show that 
the quasi-periodic oscillations (QPOs) detected with XMM-Newton originate from X41.4+60, the 
brightest X-ray source in M82. Correcting for the contributions of the unresolved sources, the QPO 
at a frequency of 55.8 ± 1.3 mHz on 2001 May 06 had a fractional rms amplitude of 32%, and the 
QPO at 112.9 ± 1.3 mHz on 2004 April 21 had an amplitude of 21%. The QPO frequency may 
possibly be correlated with the source flux, similar to the type C QPOs in XTE 1550—564 and 
GRS 1915+105, but at luminosities two orders of magnitude higher. X42.3+59, the second brightest 
source in M82, displayed a strikingly high flux of 1.4 x 10~ n ergs -1 cm -2 s -1 in the 2-10 keV band 
on 2001 May 6. A seven-year light curve of X42.3+59 shows extreme variability over a factor of 
1000; the source is not detected in several Chandra observations. This transient behavior suggests 
accretion from an unstable disk. If the companion star is massive, as might be expected in the young 
stellar environment, then the compact object would likely be an IMBH. 

Subject headings: black hole physics - accretion, accretion disks - X-rays: binaries - X-rays: 
galaxies - X-rays: individual (M82 X-l, CXOM82 J095550.2+694047=X41.4+60, 
CXOM82 J095551.1+694045=X42.3+59) 
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1. INTRODUCTION 

Ultraluminous X-ray sources (ULXs) are nonnuclcar, 
point-like X-ray sources in external galaxies with lu- 
minosities (assuming isotropic emission) above the Ed- 
dington limit for a 20M Q black hole (3 x 10 39 ergs 
s -1 ). Their fast variability and high luminosities in- 
dicate th ey may contain intermediate - mass black holes 
(TMB Hs: IColbert fc Mushotzkvl 119991 : iMakishima et all 
I2000L IKaaret et all 1200 111 However, if t he emission is 
beamed ([King et al.ll200ll: iKording et all [20021), exceeds 
the Eddington limit (IWatarai et al. 120011 : iBegelmanl 
120021 ). or both ([Poutanen et all 120071 ). IMBHs are not 
required. 

iPtak fc Griffiths! (l999) 



Using ASCA, IPtak fc Griffiths! (| 1999ft and 
iMatsumoto fc Tsurul ([19991) found that the bright 



X-ray emission from the central region of the starburst 
galaxy M82 was highly variable on time scales of hours 
to days, suggesting it arises from a compact object. 
The emission from the core of M82 was resolved with 
Chandra into several point sources, among which an 
extremely bright ULX CXOM82 J095550. 2+694047 
was identified at a lum i nosity of 10 40 — 10 41 ergs s _1 
([Matsumoto et all 120011 : IKaaret et al.ll2001ft . Following 
the convention of naming sources in M82 by their offset 
from a = 09 h 51 m 00 s , S = +69°54'00" (B1950), we 
refer to this source as X41.4+60 and use the same 
convention in referring to other X-ray and radio sources 
in M82. The source is highly variable, indicating it 
is not a supernova or remnant, and not coincident 
with the dynamical center of M82, indicating it is not 
an active galactic nucleus (AGN). If the radiation is 

1 Email: hua-feng@uiowa.edu 



isotropic and Eddington-limited, the source would be an 
IMBH candidate with a mass in excess of 500M Q . The 
luminosity of X41.4+60 is too high to be explained by 
mech anically beamed emissi on from a stellar mass black 
hole ([King fc Penned [2005) . An X-ray flare detected 
using the Rossi X-ray Timing Explorer (RXTE) in 
February 2005 showed no strong radio emission in 
contemporaneous Very Large Array ( VLA) observations 
which ruled out the possibility that t he source is 
relativistically beamed ([Kaaret et al.l[200^ ). 

XMM cannot resolve the multiple point sources de- 
tected with Chandra in the central region of M82, 
but provides superior collecti on area allowing sensitive 
studie s of timing properties. IStrohmaver fc M ushotzkv 
( 2003) analyzed a 27 ks XMM observation of the source 
obtained in 2001 and discovered quasi-periodic oscilla- 
tions (QPOs) at a frequency of 54 mHz in the energy 
range 2-10 keV. The QPOs were confirmed in a 103 ks 
XMM observat ion in 2004, however, th e frequency shifted 
to 11 4 mHz (jDewangan et all 120061 : IMucciarelli et all 
[2006ft . The power spectral density (PSD) from the sec- 
ond obse rvation also show e d a l ow frequency break at 
34 mHz. IMucciarelli et all (|2006l) found that the QPOs 
changed frequency from 107 mHz to 120 mHz during the 
second observation, and discovered a plausible harmonic 
ratio of 1:2:3 for QPO frequencies with more data from 
RXTE. Due to the limited angular resolution of XMM 
(and RXTE), all of these analyzes treated the emission 
as arising from a single point source, called M82 X-l. 

However, higher angular resolution observations made 
with Chandra show that M82 contains two ultralu- 
minous X-ray sources, CXOM82 J095550.2+694047 
(=X41.4+60) and CXOM82 J095551. 1+694045 
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(=X42.3+59), and several dimmer sources which 
contribute to the emission ascribed to M82 X-l. The 
second ULX, X42.3+59, has shown pronounced vari- 
ability, appearing as the second brightest source in M82 
on 1999 Oct 28, but being undetected on 2000 Jan 20 
(jMatsumoto et al.ll200H ; lKaaret et al.ll2001h . A Chandra 
observation in February 2005 revealed X42.3+59 in a 
bright state with strong variability, but the sensitivity 
was not adequate to confirm or reject the presence of 
QPOs. 

It is of interest to determine which of the two ULXs 
is the source of the QPO detected with XMM. We revis- 
ited the two XMM observations of M82 (g2]), corrected 
their astrometry using a Chandra observation fi j2.1[) . re- 
solved the source count rates of X41.4+60 and X42.3+59 
from surface brightness fits (£ j2.2j) . and compared them 
with timing analysis done with various source regions 
to determine the QPO origin ( 32.4|l . A seven-year light 
curve of X42.3+59 using all available XMM and Chan- 
dra archival data is presented in £ 12.51 The results are 
discussed in £J3] 

2. OBSERVATIONS AND DATA ANALYSIS 

XMM-Newton observations of M82 made on 2001 
May 06 (ObsID 0112290201) and 2004 April 21 (Ob- 
sID 0206080101) were used for image and timing analy- 
sis. The Chandra observation with the High Resolution 
Camera (HRC) made on 1999 Oct 28 (ObsID 1411-1) 
was used to correct the astrometry of XMM images. All 
Chandra and XMM observations were used to produce 
the light curve of X42.3+59. We used SAS 7.0.0 with 
up to date calibration files for XMM data reduction and 
CIAO 3.3.0.1 with CALDB 3.2.2 for Chandra data re- 
duction. 

We corrected the XMM astrometry by matching point 
sources on the Chandra and XMM images. Point sources 
in XMM images were found using edetect_chain in the 
2-10 keV band. Point sources within 18" of X41.4+60 
were not used because the diffuse emission and the two 
unresolved bright ULXs could bias the XMM source de- 
tection. Positions for Chandra sources, shown in Table_EJ 
were obtained using wavdetect, and then corrected by 
shifting X41.4+60 to the position R.A. = 09 h 55 m 50!17, 
Decl. = +69 o 40'46'.'7 (jKaaret et al.ll2001h . We note that 
none of our results depend on the absolut e astrometry. 
We a dopt the position of X41.4+60 from (jKaaret et al.l 
120011 ) as a reference in order to keep the source names 
consistent with those previously used in the literature. 
We found 8 matched point sources between the MOS1 
and HRC image for the first observation and 7 for the 
second observation. All these sources were used to align 
the XMM MOS1 images to the HRC image. 

2.1. Astrometry 

2.2. Surface Brightness Fits 

We use only XMMMOS1 images for the surface bright- 
ness fits. MOS2 images are not used because the MOS2 
point spread function (PSF) is not axisymmetric within 
a few arcseconds of the core. PN images are not suit- 
able because its pixel size of 4" slightly undersamplcs 
the core and is close to the angular distance between the 
two ULXs. M82 X-l is a bright X-ray source, so the 
MOS1 images with 3 x 10 4 source photons in 2-10 keV 



TABLE 1 

Chandra SOURCE POSITIONS USED IN THE SURFACE BRIGHTNESS 
FITS 



source # 


R.A. 


Dccl. 


name 




(J2000) 


(J2000) 




(1) 


(2) 


(3) 


(4) 


1 


09 55 54.68 


+69 41 01.1 




2 


09 55 52.31 


+69 40 54.1 




3 


09 55 51.48 


+69 40 36.0 




5 


09 55 51.05 


+69 40 45.3 


X42.3+59 


7 


09 55 50.17 


+69 40 46.7 


X41.4+60 


9 


09 55 46.61 


+69 40 41.1 





Note. — The first column indicates the source index 
in IMat sumoto ct al. (2001). All positions are obtained from 
wavdetect and then shifted by matching source 7 to the position 
R.A. = 09 h 55 m 50U7, Decl. = +69°40'46y7 JKaaret et al.ll200lT) . 



band for the first observation and 8 x 10 4 for the second 
provide adequate statistics. 

MOS1 images for the 2-10 keV band were created in 
a 250" x 250" region around X41.4+60 and sampled on 
a 1" x 1" grid with events screened by FLAG equal to 
#XMMEA_EM and PATTERN between and 12. The 
effective exposure for the MOS1 image is 30.0 ks for the 
first observation and 101.2 ks for the second. The PSF 
for MOS1 can be well described by a King function and 
the King parameters are stored in the latest calibration 
file (e.g., XRT1_XPSF.0007.CCF for MOS1 by now). We 
used the command calview with an "EXTENDED" ac- 
curacy level to produce the on-axis King model PSF at 
2 keV for MOS1. The King model parameters vary lit- 
tle with energy 2 . Therefore, a monochromatic PSF at 
2 keV provides an adequate description of the photon 
distribution for a 2-10 keV point sour ce image. 

Six of the nine sources listed in IMatsumoto et al.l 
(2001) were used in the surface brightness fit s. However, 
the source names in IMatsumoto et all ()2001l ) do not re- 
tain the sub-arcsecond accuracy that is necessary for the 
fit. We thus used the positions obtained from wavdetect 
as mentioned in £ 12.11 and listed in Table [TJ 

The 2-D surface brightness fits of the XMM MOS1 im- 
ages were performed with the CIAO program Sherpa. 
The PSF is loaded with a size of 128" x 128" around the 
core. We use 5-functions to model the point sources and 
a constant to model the background. Point sources are 
fixed at positions quoted in Table [1] and allowed to shift 
together. 

For the first XMM observation, only two point sources, 
X41.4+60 and X42.3+59, were considered in the fit. The 
model had five free parameters: the pattern shifts AX 
and AY, source amplitudes for the two ULXs (X41.4+60 
and X42.3+59), and a background amplitude taken as 
constant across the whole field. We applied the CHI 
GEHRELS statistics in Sherpa because many pixels near 
the edges of the field contained few photons. The best-fit 
source position is 0"8 away from the Chandra position, 
which implies that the astrometric correction was suc- 
cessful. The resolved count rate is 0.246 counts s^ 1 for 
X41.4+60 and 0.395 counts s' 1 for X42.3+59. The back- 
ground level is 1.6 x 10~ 6 counts s _1 arcsec -2 . In a 18"- 
radius circle, the model and data result in \ 2 = 948 with 

2 see §3 and Figure 2 in XMM-Newton current calibration file 
release notes 167. 
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TABLE 2 

Resolved 2-10 keV count hates of X41.4+60 and X42.3+59 

FROM SURFACE BRIGHTNESS FITS WITH XMM-Newton MOS1 
IMAGES 



# 


ObsID 


X41.4+60 rate 


X42.3+59 rate 


x 2 /dof 






(counts s" 1 ) 


(counts s" 1 ) 




(1) 


(2) 


(3) 


(4) 


(5) 


1 


0112290201 


0.246±0.007 


0.395±0.008 


948/1009 


2 


0206080101 


0.305±0.003 


0.064±0.003 


1453/1005 



Note. - Col. (1): Observation number. Col. (2): XMM- 
Newton ObsID. Col. (3): Best-fit source count rate of X41.4+60. 
Col. (4): Best-fit source count rate of X42.3+59. Col. (5): \ 2 and 
degrees of freedom in an 18"-radius region around X41.4+60. All 
errors are at la level. There is a systematic uncertainty on the 
count rate of X42.3+59 of 5% in the first observation and 30% in 
the second due to very nearby sources, see the text, which is not 
included in the purely statistical errors quoted in the Table. 

1009 degrees of freedom (dof), indicating an adequate fit. 



5 




-10 -5 5 



a (arcsec) 




-10 -5 5 

a (orcsec) 



Fig. 1. — 2-10 keV XMM MOS1 contour maps of the surface 
brightness in the central M82 region for the first observation (top) 
and the second observation (bottom). Coordinates are relative to 
the Chandra position of X41.4+60 (plus at the origin). The other 
plus off the origin is the Chandra position of X42.3+59. The XMM 
images have been shifted according to the best-fit positions from 
the surface brightness profile fitting. The levels are 1.5, 2.0, 2.5, 
3.0 and 3.5 (10 -3 counts s _1 arcsec -2 ) for the first observation 
and 0.7, 1.0, 1.5, 2.0, 2.5 and 3.0 (10~ 3 counts s" 1 arcsec -2 ) for 
the second. 

The situation for the second observation is slightly 
more complex, because the long exposure causes the dif- 
fuse emission and other dim sources in the central re- 
gion to become significant. If we only model X41.4+60 
and X42.3+59 plus a constant background, the x 2 /dof 



in the central 18" region is 2230/1009. We thus added 
another 4 point sources to th e mod el, i.e., sources 1, 2, 
3 and 9 in iMatsumoto et al.1 (|2001l ). Sources 4, 6 and 
8 are not added because source 8 was not bright in the 
second observation and sources 4 and 6 are too close to 
X42.3+59 (source 5). We then performed the fit with 
a model including 6 point sources and a constant back- 
ground with 9 free parameters. The 6 point sources were 
still set to move together as a whole pattern and their 
best-fit position is only 0'.'3 away from the Chandra po- 
sition. In this observation, the count rate of X41.4+60 
increased to 0.305 counts s" 1 while the X42.3+59 rate 
decreased to 0.064 counts s . The background level 
is 3.3 x 10" 6 counts s _1 arcsec -2 . Modeling the ex- 
tra 4 sources improved the fit significantly, resulting in 
X 2 /dof = 1453/1005. We checked the residual map and 
found most residuals came from structure in the diffuse 
emission region which is not included in our model. Even 
though the fit is not formally adequate, we believe the 
resolved the count rates are reliable because the emis- 
sion from the ULXs is much stronger than the diffuse 
emission. 

After checking all Chandra ACIS (Advanced CCD 
Imaging Spectrometer) observations, we found that the 
count rate in the 2-8 keV band of source 4 varied by a 
factor of 2, while source 6 appeared constant. We esti- 
mate that sources 4 and 6 contribute a systematic error 
to the count rate estimated for X42.3+59 of 5% in the 
first observation and 30% in the second. 

Contour maps of the surface brightness (corrected with 
exposure map and mask) for the two observations are 
shown in Figure[T] From the contour maps we can clearly 
see that the count rate of X41.4+60 is slightly lower than 
X42.3+59 in the first observation and much higher than 
X42.3+59 in the second. This is consistent with the rates 
from the surface brightness model fits. 

2.3. Flux from X41.4+60 

Because XMM is unable to resolve X41.4+60 from 
other point sources and the diffuse emission in M82, one 
needs great caution to make energy spectral analysis with 
XMM data. So far the only energy spectrum of X41.4+60 
which is neither contaminat ed by X42.3+59 nor strongly 
affected by pileup is from iKaaret et al.l ()2006l ) with an 
off-axis Chandra/ ACIS observation, in which the spec- 
trum is not contaminated by other sources nor suffers 
serious and uncorrectable pile-up as found in the on- 
axis Chandra/ ACIS observations. In that observation, 
X41.4+60 presented a featureless power-law spectrum 
with a photon index of 1.67. 

We adopt a p ower-law spectral index of 1.7 
(jKaaret et al.l 120061 ) with an absorption column density 
of 3 x 10 22 cm -2 , and use PIMMS to estimate the un- 
absorbed source flux in the 2-10 keV range according 
to the resolved count rates (with a factor of 0.7, be- 
cause PIMMS accepts the counting rate in a 15" re- 
gion, which contains around 70% of a 2-10 keV PSF), 
which gives 8.5 x 10" 12 ergs" 1 cm" 2 s" 1 for the first 
observation and 1.0 x 10" 11 ergs" 1 cm" 2 s" 1 for the sec- 
ond. The corresponding luminosity is 1.3 x 10 40 ergs s _1 
and 1.7 x 10 40 ergs s _1 at a distance of 3.63 Mpc 
(|Freedman et al.lll994l ) , respectively, for the first and sec- 
ond observation. The flux is not very sensitive to the 
photon index, varying only by 10% as the photon index 
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Fig. 2. — The circle indicates a 18"-radius region around 
X41.4+60. The dashed line is perpendicular to the direction con- 
necting X41.4+60 and X42.3+59. Two regions A and B are defined 
as the two half circles divided by the dashed line, and A+B refers 
to the whole circular region. 

varies between 1.5-2.3, but a decrease of the column den- 
sity to 1 x 10 22 cm -2 will bring down the unabsorbed flux 
by 25%. 

2.4. Timing Analysis 

All PN and MOS data were used for the timing anal- 
ysis in order to increase the signal to noise ratio. We 
applied FLAG equal to #XMMEA_EP and PATTERN 
between and 4 to select PN events, and FLAG equal 
to #XMMEA_EM and PATTERN between and 12 for 
MOS events. 

We examined the timing gaps in the data and back- 
ground flares before binning events into light curves. 
CCD timing gaps are stored in the good time interval 
(GTI) extension of the events file for each CCD chip. 
GTIs to exclude background flares are generated on basis 
of the quiescent intervals in the 10-15 keV light curves. 
Continuous intervals of source light curves are created 
from combined events files with a resolution of 0.5 s in 
the common intervals of all above GTIs. However, we 
found the source was so strong that the power spectral 
density (PSD) shape did not change much whether or not 
we exclude the background flares. Therefore, we decided 
not to exclude background flares and only to exclude 
CCD timing gaps. 

Every continuous piece of the light curve longer than 
512 s was divided into segments of 512 s each. Individ- 
ual PSDs were calculated with a 1024-point fast Fourier 
transform (FFT) for each 512 s segme nt. All PSDs were 
normalized to rms (|van de Klisl [19881 and averaged to 
a final one, which corresponds to an effective exposure 
of 20.5 ks for the first observation and 50.7 ks for the 
second. The PSD was rebinned linearly at v < by a 
factor of S\ and logarithmically at v > by a factor of 
&2- We set v\y — 0.1 Hz, 5\ = 3 and 62 = 1.2 for the first 
observation and 14, = 0.16 Hz, 5\ = 5 and 82 = 1.2 for 
the second. 

Three regions defined in Figure [2] were used to extract 
events: an 18"-radius circle centered on X41.4+60 di- 
vided into two half circles by a line perpendicular to the 
line connecting X41.4+60 and X42.3+59. Region A is 
the half circle on the side without X42.3+59 and region 
B is the half circle containing X42.3+59. PSDs for re- 
gions A+B, A and B are presented in Figure [3] for the 



first observation and in Figure [4] for the second. 

We fit the PSDs for the first observation with a power- 
law plus Lorentzian model; the former component is to 
model the continuum and the latter is for the QPO. 
For the second observation, the continuum component 
is modeled by an exponentially cutoff power-law model. 
The best-fit parameters including the power-law slope a, 
cutoff frequency v m t, Lorentz centroid and full width 
at half maximum (FWHM) u>l, as well as the fractional 
rms amplitude (rms/mean) of the QPO component, are 
listed in Table [3] with la errors, respectively for each 
region and observation. We note that the QPO is not 
detected in region B for the first observation. We fixed 
the power-law index to —1.5 to model the continuum and 
evaluated a 2a upper limit of the QPO rms by integrating 
the powers in the 30-70 mHz frequency region. 

Given the resolved source count rates from the sur- 
face brightness fits, X41.4+60 contributes 41%, 67% and 
30% of 2-10 keV photons respectively for region A+B, A 
and B for the first observation, and 84%, 94% and 76% 
respectively for the second. Assuming the QPO comes 
from one of the ULXs, the required rms for a particular 
source is obtained by normalizing the measured rms to 
the photon fraction in each region. The calculated values 
are shown in columns (8) and (9) in Tabled Column (8) 
shows consistent rms values from the different regions for 
each observation. Thus, the assumption that X41.4+60 
produces the QPOs gives consistent results regardless of 
the region chosen to measure the QPO strength. Con- 
versely, column (9) shows that if X42.3+59 is assumed to 
produce the QPOs then widely differing rms values are 
found from the different regions. Thus, we conclude that 
the QPOs originate from X41.4+60. 

2.5. A Seven-Year Light Curve of X42.3+59 

All fourteen available Chandra observations (see Ta- 
ble @} and the two XMM observations were used to pro- 
duce a light curve of the 2-10 keV unabsorbed flux from 
X42.3+59 which covers a period of 7.3 years. The XMM 
observations are described above. The resolved count 
rates of X42.3+59 are 0.395 counts s _1 on MJD (modi- 
fied Julian Date) 52035.6 and 0.064 counts s" 1 on MJD 
53117.5. 

The source count rate or flux in two Chandra observa- 
tions has been reported in literature. Thes e are Chandra 
obser vation 10, see Table 01 with A CIS ([Kaaret et alJ 
2006), and observation 3 with HRC (jMatsumoto et ail 
2001 . 

For ACIS observations, we calculated the pileup frac- 
tion following the definitions in The Chandra ABC Guide 
to Pileup 3 where ft is the fraction of good events lost due 
to grade or energy migration, f e is the fraction of single 
events over all detected events, f r is the fraction of count 
rate lost, and a is the probability that a piled event is 
retained as a good grade. 

Two ACIS observations, 11 and 12, suffer only mild 
pileup with f t < 10%. X42.3+59 has the same count 
rate in these two observations. Observation 12 is about 
one day later than observation 11, and has an expo- 
sure of only 1/4 of observation 11. We assume the 
source stayed in the same state for both observations 
and use observation 11 to constrain the spectral pa- 

3 http:/ /asc. harvard.edu/ciao/download/doc/pileup_abc.ps 
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TABLE 3 

Best-fit PSD parameters for light curves extracted from different source regions in the 2-10 keV range. 



region 
(1) 


a 
(2) 


'•'cut 

(mHz) 
(3) 


1-t, 
(mHz) 

(4) 


wl X 2 /dof 
(mHz) 
(5) (6) 


rmsqpo 
(7) 


rms4i.i 

(8) 


rms42.2 
(9) 


First observation 


A+B 
A 
B 


-1.6 + 0.6 
-1.2 + 0.5 
— 1.5 fixed 




52 + 3 
55.8 ± 1.3 


29 + 9 21/25 
19 + 4 26/25 
38/29 


12.1 + 0.8% 
21.6 + 0.8% 
< 11% 


30 ± 2% 
32.2 ± 1.2% 
<37% 


20.5 ± 1.4% 
65 ± 2% 
<16% 


second observation 


A+B 
A 
B 


0.4 + 0.4 

3+ - 4 
u -°-0.3 

6 +0 ' 9 
u - D -0.6 


0.03 ± 0.02 

9+3 

A -1 
o+14 
d -3 


112.9 ± 1.3 
113 ± 2 
115 ± 3 


30 + 5 14/20 
28 + 8 11/20 
30 + 9 15/20 


17.9 ± 0.3% 
20.3 ± 1.1% 
15.6 ± 0.9% 


21.3 ± 0.4% 
21.6 ± 1.2% 
20.5 ± 1.2% 


111.9 ± 1.9% 
338 ± 18% 
65 ± 4% 



Note. — Col. (1): Regions from which the light curves are extracted; see Figure[2]for details. Col. (2): Exponent of the power-law (first 
observation) or the exponentially cutoff power-law (second observation). Col. (3): Cutoff frequency of the exponentially cutoff power-law 
(only for second observation). Col. (4): Centroid of the Lorentzian. Col. (5): Full width at half maximum of the Lorentzian. Col. (6): x 2 
and degrees of freedom. Col. (7): Fractional rms amplitude of the QPO (= the Lorentzian) Col. (8): Required rms if the QPO originates 
from X41.4+60. Col. (9): Required rms if the QPO originates from X42.3+59. All errors are at 1 a level. 
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Fig. 3. — 2-10 keV PSD for the first observation calculated 
from regions A+B, A and B. The dashed line is the best-fit model 
including a power-law component and a Lorentzian component. 
The PSD for region B is modeled only by a single power-law. 
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Fig. 4. — 2—10 keV PSD for the second observation calcu- 
lated from regions A+B, A and B. The dashed line is the best-fit 
model including a exponentially cutoff power-law component and 
a Lorentzian component. 
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TABLE 4 
Chandra Observations of M82. 



index 
(1) 


MJD 

(2) 


ObsID 

(3) 


instrument 
(4) 


exposure 
(ks) 
(5) 




51441.7 


361 


ACIS-I 


33.3 


2 


51442.0 


1302 


ACIS-I 


15.5 


3 


51479.4 


1411-1 


HRC-I 


36.3 


1 


51542.1 


378 


ACIS-I 


4.1 


5 


51563.7 


1411-2 


HRC-I 


17.8 


6 


51614.9 


379 


ACIS-I 


8.9 


7 


51671.9 


380-1 


ACIS-I 


3.9 


8 


51707.6 


380-2 


ACIS-I 


1.2 


9 


52443.9 


2933 


ACIS-S 


18.0 


10 


53406.3 


6097 


ACIS-S 


52.8 


11 


53599.5 


5644 


ACIS-S 


68.1 


12 


53600.8 


6361 


ACIS-S 


17.5 


13 


54109.7 


8189 


HRC-S 


61.3 


11 


54112.6 


8505 


HRC-S 


83.2 


Note. 
vations, 


— ObsID 380 and 1411 consist of two individual obser- 
respectively. The MJD indicates the mean time of each 



observation. 



The source was not detected in five observations: three 
with HRC (5, 13, and 14) and two with ACIS (1 and 
2). For the HRC observations, the upper limit of the 
count rate is calculated from a 1" source region. For the 
ACIS observations, it is adopted from the brightest pixel 
around the source region. We note that there is strong 
diffuse emission around X42.3+59, thus our upper limits 
are conservative and likely an overestimate. 

All count rates (with a factor of 0.7 for XMM obser- 
vations) are converted to 2-10 keV unabsorbed flux with 
PIMMS given Nn = 3 x 10 22 cm~ 2 and a power-law pho- 
ton index T — 1.5. The flux estimates from the XMM 
data include the systematic uncertainty caused by the 
very nearby sources 4 and 6 as mentioned above. The 
complete light curve with all the available XMM and 
Chandra data is presented in Figure O It is clear that 
X42.3+59 is highly variable. The source exhibited an 
unabsorbed flux of 1.4 x 10~ n ergs -1 cm -2 s -1 and a 
corresponding luminosity of 2.2 x 10 40 ergs s _1 in the 
first XMM observation, reaching its brightest state that 
has ever been observed. 
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1000 1500 2000 2500 3000 3500 4000 4500 
MJD - 50000 

Fig. 5. — A seven-year light curve of X42.3+59 observed with 
XMM/MOS, Chandra/ ACIS or Chandra/URC. The left axis in- 
dicates the 2—10 keV unabsorbed flux, and the right axis is the 
corresponding luminosity at a distance of 3.63 Mpc. 

rameters. The energy spectrum is created from a 1" 
source region binned for a minimum counts of 100 per 
bin. We subtracted a background estimated using a 
nearby diffuse emission region free of point sources and 
then fitted with a pileup corrected, absorbed power-law 
model in the energy range of 0.3-8 keV using Sherpa. 
Best-fit parameters include an absorption column density 
Nu = (3.45 ± 0.11) x 10 22 cm" 2 and a power-law index 
T = 1.49±0.06, with X 2 = 129 for 101 degrees of freedom. 
The pileup model indicates a pileup fraction f e = 1.2% 
with a = 0.28, which is well consistent with the estimate 
from the count rate. The best-fit model predicts an un- 
absorbed source flux of 6.0 x 10~ 12 ergs -1 cm -2 s _1 in 
2-10 keV. 

Five ACIS observations (4, 6, 7, 8, and 9) suffer from 
severe pileup. These observations have a pileup fraction 
ft > 20% and, thus, are not suitable for spectral analysis. 
The count rates are measured from a 2 x 2 pixel array in 
0.3-8 keV and then corrected by f r , which is estimated 
from a 3 x 3 pixel island in the full band assuming a = 0.5. 
The background level, measured from observation 11, is 
only about 2% of the source flux for these observations 
and thus negligible. 



3. DISCUSSION 
3.1. Origin of the QPO 

We present solid evidence that the QPOs from M82 
originate from the brightest source X41.4+60, rather 
than the second brightest source X42.3+59. The angu- 
lar resolution of XMM is not adequate to cleanly resolve 
the two sources, but is adequate, with surface brightness 
fitting using the known source positions, to determine 
the count rate from each source. Using these count rates 
and the QPO amplitudes calculated from three differ- 
ent extraction regions, we find that the assumption that 
X41.4+60 produces the QPOs give consistent results. 
Conversely, if X42.3+59 is taken as the source of the 
QPOs, then the QPO rms amplitude values calculated 
from different extraction regions differ at a significance 
level of at least 20er. Therefore, we conclude that the 
detected QPOs originate from X41.4+60 at a frequency 
of 55.8 mHz with an rms of 32% in the 2-10 keV band 
in the first observation, and at a frequency of 112.9 mHz 
with an rms of 21% in the second observation. 

The origin of the broken p ower-law feature i n the PSD 
continuum, first reporte d bv lDewangan et ail (|2006f ) and 
iMucciarelli et al. (2006) in the second XMM observation, 
has not yet been determined. The same feature is con- 
firmed in our analysis (although we fit this noise compo- 
nent with an exponentially cutoff power-law). The inte- 
grated powers for the continuum in the 6-80 mHz range 
are 14±4%, 16±6% and 11 ± 5%, respectively for region 
A+B, A and B. Assuming only one of the sources pro- 
duces the continuum noise, we require an rms amplitude 
of 17 ± 5%, 17 ± 6%, and 14 ± 7% from X41.4+60, or 
87 ± 25%, 267 ± 100%, and 46 ± 21% from X42.3+59, 
to produce the detected continuum noise in region A+B, 
A and B, respectively. Because the errors are large, we 
cannot completely rule out X42.3+59 as the origin of the 
continuum noise. However, it appears that X41.4+60 is 
more likely to be the origin of the continuum power. 

3.2. X41.4+60 

It is of interest to use the detected QPO frequencies 
to attempt to constrain the mass of the compact object 
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Fig. 6. — 2-10 keV unabsorbed power-law luminosity versus 
QPO frequency for XTE J1550-564 (circle), GRS 1915+105 (tri- 
angle), X41.4+60 (square) and NGC 5408 X-l (star). 

in X41.4+60. Provided the QPO frequency is limited by 
the Keplerian frequency at the innermost stable circu- 
lar orbit around a Schwarzschild black hole, the compact 
object mass of X41.4+60 is constrained to be less than 
3 x 1O 4 M0 according to the maximu m QPO frequency o f 
190 mHz that has been reported in lKaaret et al.l (|2006f) . 
This clearly rules out that X41.4+60 could be a super- 
massive black hole. 

The centroid of low frequency QPOs in Galactic black 
holes varies in a large range , e.g. , 10 -2 -10 Hz for 
GRS 1915+105 ijMorgan et al]ll997h . Therefore, it is 
difficult to derive the compact object mass by simply 
scaling the QPO frequency to mass between ULXs and 
Galactic black holes. One must, instead, determine if the 
QPO frequency correlates with another observable, such 
as the luminosity or spectral state, which would enable 
one to correct for the intrinsic frequency variations. 

The X41.4+60 QPO in the second XMM observa- 
tion is analogous to t he "type C" QPO as defined by 
iRemillard etail ((2002): appearing with a flat-top PSD 
and high amplitude. The QPO in the first XMM ob- 
servation is also of type C that has a high ampli- 
tude. A positive correlation between the QPO fre- 
quency and disk/power-law flux has been found for type 
C QPOs in some stellar mass black hol e binaries like 
XTE J1550-564 and GRS 1915+105 (iSobczakeTall 
I2000al iReig et al.l 120001: iRemillard fc McClintockl l2006fh 
We note that GRO J1655— 40 displayed a negative cor- 
relat ion between the QPO frequency and the power-law 
flux (|Sobczak et al.ll2000al) . However, those QPOs have 
relativel y small amplitudes an d Q factors are likely not of 
type C |Remillard et al.lfl999h . We obtained timing and 
spect ral data for two Galac tic sources: XTE J550— 564 
fromlSobczak et al.lj|2000al |b1) and GRS 1915+105 from 
IVignarca et~aH (|2003l ). For XTE J550-564, we consider 
only QPOs with amplitudes ^13% to include only type 
C QPOs. The QPO frequency versus 2-10 keV luminos- 



ity is presented in Figure [6] For comparison, we plotted 
the data for X41.4+60. We include an uncertainty of 
25% on the XMM luminosities for X41.4+60 to account 
for uncertainties in the conversion from XMM-Newton 
count rates to luminosities due to the lack of knowledge 
of the spectrum. The two observations weakly suggest a 
similar positive correlation between the QPO frequency 
and the source flux (power-law flux). However, due to 
the significant uncertainty in the luminosities extracted 
from the XMM-Newton data and the fact that there are 
only two data points, any conclusions must be regarded 
as tentative. We also plotted one data point for the ULX 
NGC 5408 X-l, which is the second ULX found to pro- 
duce QPOs, in this case near 20 mHz (|Strohmaver et all 
[20071) . The QPOs detected from NGC 5408 X-l are also 
of type C. 

GRS 1915+105 is slightly more massive than 
XTE J1550— 564, and its luminosity at a given QPO fre- 
quency is slightly higher than XTE J1550-564. If the 
frequency versus luminosity pattern does, indeed, scale 
with the mass of the compact object, then the data in 
Figure [6] would indicate that X41.4+60 and NGC 5408 
X-l are IMBHs. Though the sample of QPO detections 
for the ULXs is not large enough to derive the black hole 
mass accurately, we can see in Figure [6] that the pattern 
for X41.4+60 is shifted by a factor of 100-500 in luminos- 
ity relative to that of the Galactic black holes, indicating 
a black hole mass of around 10 3 M Q . More observations 
of both ULXs are required to check if the QPO frequency 
robustly follows a correlation with luminosity similar to 
that seen from Galactic stellar mass black hole X-ray 
binaries. 

Correlations between the QPO and the power-law 
photon ind ex have been observed in several Galacti c 
black holes (ISobczak et alll2000at IVignarca et al.ll2003l) . 
iTitarchuk fc Fioritol (|2004l ) have suggested that the QPO 
frequency scales inversely with the mass of the compact 
object at fixed power-law photon index. Unfortunately, 
it is difficult to accurately measure the photon index for 
X41.4+60 or X42.3+59 with XMM due to large overlap 
between their PSFs. This prevents us from making a 
quantitative comparison with the relations seen for stel- 
lar mass black hole X-ray binaries. It would be of inter- 
est to conduct simultaneous XMM and Chandra observa- 
tions in order to derive timing information from XMM si- 
multaneously with spectral information from Chandra in 
order to study the QPO frequency versus spectral index 
correlation for X41.4+60. A combined spectral/timing 
analysis of the behavior of X41.4+60 would likely pro- 
vide the most robust estimate of compact object mass. 
The same observations could also greatly improve the 
uncertainty in luminosity and, thus, provide a test of 
the putative QPO frequency versus luminosity correla- 
tion discussed above. 

The low frequency break (from slope to —1) in 
the PSD is tightly correlated w ith the QPO frequency 
(fWijnands fc van der Kl"islll999[ K We note the 34 mHz 
break frequency found in M82 X-l (plausibly from 
X41.4+60) and the 113 mHz QPO frequency are con- 
sistent with this correlation. The break frequency in the 
PSD has been thought to be correlated with the mass of 
the compact object fe.g.. iMarkowitz et al.l [2003h . Com- 
paring to Cygnus X-l, which has a low frequency break 
varying between 0.02-0.4 Hz (|Belloni fc hasingeri fl990L 
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iNowak et allll999D . the compact object mass of X41.4+60 
would be constrained to be betwee n 6-120 M^, assum - 
ing a mass of 10M Q for Cygnus X-l dHerrero et alJfl995h . 
We note that the estimated mass is lower than that in- 
ferred from isotropic emission under the Eddington limit, 
or the QPO frequency-luminosity correlation discussed 
above. Again, we emphasize the importance of obtaining 
multiple observations simultaneously with Chandra and 
XMM-Newton to study the correlations between spec- 
tral and timing properties which could lead to a robust 
estimate of the compact object mass. 

3.3. X42.3+59 

The seven-year light curve of X42.3+59, see Figure El 
shows that the ULX is highly variable. The peak lumi- 
nosity implies the compact object mass is at least 200M Q 
assuming isotropic emission. The ratio of the maximum 
to minimum observed flux is at least 1000. 

The huge variability of X42.3+59 could be a natural 
consequence of a relativistic jet aimed along our line of 
sight - a so called "microblazar" . X42.3+59 is coinciden t 
with the radio source 42.21+59.0 (|Muxlow et al.lll994D . 
However, the radio counterpart presents a thermal spec- 
trum instead of a synchrotro n spectrum (jMuxlow et alJ 
11994 iMcDonald et alJ l2002ft shows little, if any, vari- 



ability on long time scales rtKronberg et al.ll200C), an d 
has a spatial extent of 4.9 pc ( McDonald et alJ "2002). 
These results indicate that X42.3+59 is not relativistic 
je t emission. 

iKalogera et al.1 (|2004h suggested an observational test 
to distinguish whether ULXs are stellar mass objects 
with beamed emission or IMBHs. The former are likely 
to arise as thermal time-scale mass transfer binaries, in 
which the high mass accretion rates lead to thick accre- 
tion disks causing geometrically beamed emission. Ther- 
mal time-scale mass transfer generally produces stable 
disks and persistent X-ray emission. X-ray binaries are 
transient when the disk temperature at outer edge is be- 
low the hydrogen ionization temperature. This requires 
an average ma ss transfer rate be low a critical value. As 
first shown bv King et al.l (|1996h . this implies that the 
black hole mass must be above a minimum set by the 
companion star mass and the orbital period. Transient 
behavior with a massive companion star likely requires 
anIMBH. 

Unfortunately, the nature of the X42.3+59 binary is 
unknown. Identification of the spectral type of the com- 
panion star and measurement of the orbital period, to- 



gether with the knowledge that the system is an X-ray 
transient, could be used to place direct constraints on 
the mass of the compact object. X42.3+59 is thought 
to be associated with an Hn region (|Muxlow et al.|[l"994h 
th at is full of young, ho t stars. Simulations (Fig. 1 and 2 
in IKalogera et al. 2004) show that if the companion was 
initially a 10-20M© star, then transient behavior, such 
as that observed, likely requires an IMBH. 

From the light curve in Figure El it seems that 
X42.3+59 has just completed an outburst that lasted 
about 6 years. This outburst duration is longer than 
that typically seen from Galactic soft X-ray transients, 
but could be similar to the current outburst of the tran- 
sient GRS 1915+105. Given an accretion rate (M) and 
donor mass (Afg), th en according to equation (2) of 
IKalogera efaLl (|2004h . an IMBH (> 2OM ) is required 
for a transient system if the binary period is less than a 
threshold value 



Pu 



2.4 yr 



AI 



10- 4 M G yr- 



M 2 



10 M, 



+3 



(1) 



The average accretion rate M can be estimated from 
the outburst luminosity L ont assuming an outburst duty 
cycle d and an accretion power efficiency e as M = 
dL out / (ec 2 ). We can, thus, rewrite the equation above 
as 



6.8 d x 



d L fl 



e 10 40 ergs s _1 



M 2 
1 Mr. 



(2) 



From Figure El the outburst luminosity of X42.3+59 is 
about 10 40 ergs s _1 . Assuming d — 0.1 and e = 0.1, then 
M ~ 1.8 x 10~ 7 Af Q yr -1 . If the companion mass Mi > 
1M©, then the threshold period for transient behavior 
for a black hole mass larger than 20M Q is 6.8 d. The 
threshold period depends very weakly on the companion 
mass. Measurement of the orbital period of X42.3+59 
should enable us to place a lower bound on the compact 
object mass. 
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